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^ \ The longitudinal structure function in deep inelastic scattering is one of the observables 

^ \ from which the gluon distribution can be unfolded. Consequently, this observable can be used 

^ ' to constrain the QCD dynamics at small x. In this work we compare the predictions of distinct 

^ ■ QCD models with the recent experimental results for Fl{x, Q^) at small x and low Q^ obtained 

^ ! by the HI collaboration. We focus mainly on the color dipole approach, selecting those models 

^ \ which include saturation effects. Such models are suitable at this kinematical region and also 

O \ resum a wide class of higher twist contributions to the observables. Therefore, we investigate 

^ ' the influence of these corrections to F^ in the present region of interest. 
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H 1 Introduction 

c^ , 

The small x regime in deep inelastic scattering (DIS) is one of the frontiers of Quantum 
Chromo dynamics (QCD). This represents the challenge of studying the interface between 
perturbative and nonperturbative QCD, with the characteristic feature that the transition 
is taken in a kinematical region where the strong coupling constant Ug is small. That re- 
gion has been explored by the electron-proton (ep) scattering at HERA, which has shown a 
striking rise of the proton structure function F2{x,Q'^) for values x < 10~^. This behavior 
implies that the cross section increases faster that logarithmically with the energy, violat- 
ing the Froissart bound. Therefore, new dynamical non-linear QCD effects associated to the 
unitarity corrections are expected to slow down its further growth [1, 2, 3, 4, 5, 6, 7]. The 
search of signatures for these effects has been an active subject of research in the last years 
[8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24]. 



In particular, it has been observed that the HERA data at small x and low Q^ can be 
successfully described with the help of saturation models [10, 11, 12, 21]. Moreover, the exper- 
imental results for the total cross section [14] and also for the inclusive charm production [20] 
present the property of geometric scaling, which is one of the main characteristics of the high 
density QCD approaches (For a recent review see e. g. Ref. [25]). The saturation (non-linear 
QCD) approaches are characterized by a typical scale, denoted the saturation scale Qs{x), 
which is energy dependent, and marks the transition between the linear (leading twist) pertur- 
bative QCD regime and saturation domain. As current phenomenological saturation models 
have indicated that for the HERA domain the saturation scale is smaller than 2 GeV^, we 
expect that the signatures of the saturation effects becomes more evident in the region of small 
X and very low Q^. Furthermore, some of these approaches contain information of all orders 
in l/Q"^, namely they resum higher twist contributions [26, 27]. These corrections should be 
important at the low Q^ region, where the leading twist (DGLAP) approaches would be in 
the limit of their aplicability. Therefore, saturation models are quite suitable for the present 
phenomenological study of the longitudinal structure function at low Q^. 

In this work we analyze the behavior of the longitudinal structure function Fl{x, Q"^) in this 
kinematical regime. One considers several QCD theoretical approaches, focusing mainly on the 
saturation models. The predictions are compared with the recent (preliminary) Fl experimental 
results, as determined from the 1999 minimum bias and the 2000 shifted vertex HI data [28]. 
A comment related to these results is in order here. The experimental determination of Fl 
is difficult since it usually requires cross sections measurements at different values of center 
of mass energy, implying a change of beam energies. An alternative possibility is to apply 
the radiation of a hard photon by the incoming electron. Such hard radiation results into an 
effective reduction of the center of mass energy. Several studies on the use of such events to 
measure F^ have been carried out [29]. With these measurements, which in principle could be 
performed in the near future, it may be possible to explore the structure of Fl{x, Q^) in the low 
x range. However, currently, to obtain the F^ data, the HI Collaboration has parametrized the 
structure function F2 taken only data for y < 0.35, where the contribution of Fl is small. This 
parameterization was evolved in Q^ according to the DGLAP evolution equations, providing 
predictions for F2 in the high y region which allowed, by subtraction of the contribution of F2 to 
the cross section, the determination of the longitudinal structure function (For more detailed 
discussions see Refs. [30, 28]). Therefore, the Fl data only are obtained after the use of a 
procedure in the measurements of the total cross section. 

The paper is organized as follows. In the next section, we briefly review the theoretical 
description of Fl in the linear DGLAP (leading twist) approximation and summarize the main 
expressions considering the color dipole approach. For the latter, we introduce two representa- 
tive saturation models which have their phenomenological parameters well constrained from the 
current small x experimental data. The comparison of the numerical results, further discussions 
and conclusions are presented in the last section. 



2 Theoretical description of Fi at small x 

The longitudinal structure function F^ corresponds to the interaction of the longitudinally 
polarized virtual photon in the one-photon-exchange mechanism of lepton-nucleon scattering. 
It is a very interesting dynamical quantity since, at least at low x, its dominant contribution 
comes from gluons. While in the naive parton model this structure function vanishes, at leading- 
order in asiQ"^) it acquires a leading twist contribution. At small x this contribution is driven 
by the gluon through the g ^ qq transition and, in fact, F^ can be used as a very useful quantity 
for a direct measurement of the gluon distribution in a nucleon. One can write the longitudinal 
structure function Fl in terms of the cross section for the absorption of longitudinally polarized 
photons as 

FLix, Q') = ^iin^ at'ix, Q') ^ -^ af'ix, Q') (1) 

at small x. Therefore, accurate measurements of Fl at low x and/or Q^ would be helpful to con- 
strain the physics in that kinematical region. In particular, we expect that this observable may 
discriminate between the leading twist predictions, which consider the collinear factorization 
and parton distributions determined from global fits, and the predictions from the saturation 
models which resum a class of higher twist contributions at small x. 

While the longitudinal structure function is (at least theoretically) fairly well understood at 
high Q^ very little (if anything) is known about its possible extrapolation towards the region 
of low Q^ and small x [31] (For recent discussions see e. g. [32, 33]). Theoretically, we have 
that in the limit Q^ — >■ the structure function Fl has to vanish as Q"^. It reflects the simple 
physical fact that the total cross section ctl ~ Fl/Q^ describing the interaction of longitudinally 
polarized virtual photons has to vanish in the real photoproduction limit. On the other hand, 
the leading twist DGLAP MRST [34] and CTEQ [35] global fits require the gluon distribution 
to be valencelike or negative at small x and low Q^ in order to describe the experimental 
data, leading to Fl being negative at the smallest x — Q^. At that region, a comparison of 
the predictions at LO, NLO and NNLO using MRST partons shown a poor description of the 
experimental results [36]. However, the description is improved if a ln(l/a;) resummation is 
considered [37]. Here, it is important to emphasize that higher-twist contributions are not 
considered in these analyzes. Nevertheless, in the global fit of the existing light-targets DIS 
data at LO, NLO and NNLO QCD approximations, Alekhin [38] has estimated the high-twist 
contributions to the structure functions. It was verified that these terms do not vanish up to 
NNLO (See Fig. 12 in that reference) and gives important contributions at both small and large 
X regions. Although the expectation that higher twist plays an important role at very large x is 
not new, the contribution of these terms in the small x region has been a subject of discussion 
only in the last years. For instance, in Ref. [39] a simple parameterization of the higher twist 
contribution to the F2 structure function have been used, and it was found that for x < 0.5 
the resulting correction is small and negative but beyond 0.6 large and positive. In particular, 
the higher twist contributions for x < 0.01 are very small at HERA low x domain. In other 
words, the experimental results for the F2 structure function in principle can be described by 
a leading twist approximation. However, this feature can also be explained as being due to 



the almost complete cancellation of the twist-4 corrections to the transverse and longitudinal 
structure functions [26]. Therefore, only a direct analyzes of F^ could discriminate between 
leading twist and higher-twist resummations. 

In order to address these issues, in what follows we present representative theoretical ap- 
proaches taking into account the usual DGLAP leading twist approximation and the twist 
resummation rendered by the saturation models. At leading order, twist-two and in the in- 
finite momentum frame the longitudinal structure function can be expressed in terms of the 
Altarelli-Martinelli equation [40] 
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which shows the dependence of Fl on the strong constant coupling and on the gluon density. 
At small X, the second term is the dominant one since it is driven by the gluon distribution. 
Consequently, Eq. (2) can be reasonably approximated by F^ ^ 0.3 ^ x g (2.5 x, Q'^) [41]. This 
relation demonstrates the close relation between the longitudinal structure function and the 
gluon distribution. In our further numerical calculations using the Altarelli-Martinelli equation, 
one considers as input the MRST2001(LO) [34] and GRV98LO [42] parton distributions. We use 
the GRV98 parameterization in order to compare the coUinear approach with the experimental 
results for Q^ < 1 GeV^. 

In the proton rest frame, the DIS process can be seen as a succession in time of three 
factorisable subprocesses: i) the photon fluctuates in a quark- ant iquark pair with transverse 
separation r± ~ 1/Q long after the interaction, ii) this color dipole interacts with the pro- 
ton target, iii) the quark pair annihilates in a virtual photon. The interaction '~f*p is further 
factorized in the simple formulation [43], 
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where z is the longitudinal momentum fraction of the quark, x ~ Q'^/W^ is equivalent to the 
Bjorken variable. The photon wavefunctions '^l,t are determined from light cone perturbation 
theory and read as 
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where the auxiliary variable e^ = z{l — z) Q"^ + mj depends on the quark mass, rrif. The Kq^i 
are the McDonald functions and the summation is performed over the quark flavors. 

The dipole hadron cross section adip contains all information about the target and the strong 
interaction physics. There are several phenomenological implementations for this quantity 
[9, 10, 11, 12, 19, 21]. The main feature of these approaches is to be able to match the soft (low 
Q^) and hard (large Q^) regimes in an unified way. In the present work, we follow the quite 



successful saturation models [10, 21], which interpolates between the small and large dipole 
configurations, providing color transparency behavior, a dip ^ r^, as r ^ Qs and constant 
behavior at large dipole separations r < Qg. It is important to emphasize that in the dipole 
models at small x, both Fl and F2 are governed by adip and therefore behave similarly. In 
particular, Fl should go to zero when Q^ -^ at low x in the dipole picture since \^ l\^ oc 
Q^. The parameters of the saturation models have been obtained from phenomenological 
adjustments to small x HERA data. As a first model, we present the analytically simple 
GBW model, which resembles the main features of the Glauber-Mueller resummation. Its 
phenomenological aplication has been successful in a wide class of processes with a photon 
probe (DIS, diffractive DIS, Deeply Virtual Compton Scattering, heavy-quark production, two- 
photon physics) [10, 12, 13, 14, 15, 16, 17, 18, 20]. The parameterization for the dipole cross 
section in this model takes the eikonal-like form. 



c^dipix, r^) = (To 
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Q^(a:) = (|) GeV^ (4) 



where the parameters were obtained from a fit to the HERA data producing o"o = 23.03 (29.12) 
mb, A = 0.288 (0.277) and xq = 3.04-10^^ (0.41 ■ 10"^) for a 3-fiavor (4-fiavor) analysis [10]. An 
additional parameter is the effective light quark mass, rrif = 0.14 GeV, which plays the role of 
a regulator for the photoproduction {Q^ = 0) cross section. 

An important aspect of the saturation models is that they resum a class of higher twist 
contributions which should be non-negligible in the low Q^ regime [26, 27]. Consequently, in 
this kinematical region we may expect a discrimination between the twist two calculations, 
usually considered in the global fits of the experimental data, and the saturation models. Some 
hints of the differences between these models has been presented in Refs. [26, 27]. In particular, 
the twist expansion of the GBW model has been calculated, with the different twist terms in 
the massless limit given by: 

(4) _ Y^ 2'^em I y^ r^2 I ^ I f'i\C'^ C2 1 
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where ^ = ^ is the scaling variable which appears in the geometric scaling property of the 
inclusive cross section and iIj{x) is the digamma function. The results above can be contrasted 



to the leading twist result a}^ ~ <^oJ2f^}^^^^^- 



The longitudinal twist-4 and twist-8 terms give sizeable negative corrections to the leading 
twist contribution, mainly at ^ = % ~ 1. It is expected that a precise low Q^ measurements of 
Fl at small-a: could reveal this important feature. This is investigated in the analysis presented 
here, considering the recent HI preliminary data on the longitudinal structure function. It 
should be noticed although F2 had been measured in this region with accurate precision, its 



longitudinal and transverse twist-4 contributions have opposite signs and almost the same order 
of magnitude. Hence, they approximately cancel each other and produce either a small twist-4 
correction. Therefore, for the inclusive F2 structure function the higher twist corrections are 
hidden in the mismatch between the longitudinal and transverse higher twist corrections. 

Despite the saturation model to be very successful in describing HERA data, its functional 
form is only an approximation of the theoretical non-linear QCD approaches. On the other 
hand, an analytical expression for the dipole cross section can be obtained within the BFKL 
formalism. Currently, intense theoretical studies has been performed towards an understanding 
of the BFKL approach in the border of the saturation region [44, 45]. In particular, the dipole 
cross section has been calculated in both LO and NLO BFKL approach in the geometric scaling 
region [46]. It reads as. 



CTdipK.x, r) = o-Q 
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where ctq = 27ri?p {Rp is the proton radius) is the overall normalization and the power 7sat is 
the (BFKL) saddle point in the vicinity of the saturation line Q^ = Qg^^^{x). In addition, the 
anomalous dimension is defined as 7 = 1 — 7sat- As usual in the BFKL formalism, ag = Nc a.s/'^, 
j3 ^ 28^(3) and Y = ln(l/x). The quadratic diffusion factor in the exponential gives rise to 
the scaling violations. 

The dipole cross section in Eq. (8) does not include an extrapolation from the geometric 
scaling region to the saturation region. This has been recently implemented in Ref. [21], where 
the dipole amplitude Af{x, r) = adip/lnR^ was constructed to smoothly interpole between the 
limiting behaviors analytically under control: the solution of the BFKL equation for small 
dipole sizes, r <^ 1/Qsat(a;), and the Levin- Tuchin law [47] for larger ones, r 3> 1/Qsa.t{,x). 
A fit to the structure function ^2(0;, Q^) was performed in the kinematical range of interest, 
showing that it is not very sensitive to the details of the interpolation (For a comprehensive 
phenomenological analysis of the HERA results using the numerical solution of the BK equation 
see Ref. [48]). The dipole cross section was parametrized as follows. 



2 7sat + 



ln(2/rQ,^t) 



a^ix, r) = a,{ ^0 (^) ^"^' ^'"^ ^ for rQ,..{x) < 2 , (9) 

L - exp^" 1°' (6^Q=at) ^ for rQ^^tix) > 2 , 

where the expression for rQsa,t{x) > 2 (saturation region) has the correct functional form, as 
obtained either by solving the Balitsky-Kovchegov (BK) equation [5, 6], or from the theory of 
the Color Glass Condensate (CGC) [25]. Hereafter, we label the model above by CGC. The 
coefficients a and b are determined from the continuity conditions of the dipole cross section 
at rQss.t{x) = 2. The coefficients 7sat = 0.63 and k, = 9.9 are fixed from their LO BFKL 
values. In our further calculations it will be used the parameters Rp = 0.641 fm, A = 0.253, 
xq = 0.267 X 10~^ and A/q = 0.7, which give the best fit result. A large x threshold factor 
{1 — x)^ will be also considered, for sake of completeness. 

Recently, this model has also been used in phenomenological studies of the vector meson 
production [22] and the diffractive processes [23] at HERA as well as hadron production in 

6 



a 



O 
x" 




Figure 1: The results for Fl{x,Q'^) as a function of x at fixed low Q^ values. The numerical 
results for the saturation models (GBW and CGC) as well as leading twist DGLAP approach 
for three inputs for the gluon distribution (GRV98, EHKQS and MRST) are presented. Data 
from HI Gollaboration. 

nuclear collisions at RHIC [24]. Here we compare for the first time this model with the recent 
HI data for the longitudinal structure function. 



3 Results and Discussions 



Lets present the numerical results coming out the saturation models (GBW and CGC) and from 
the usual coUinear approach. In Figs. 1 and 2 the predictions of various theoretical models 
are compared with the HI experimental results for the longitudinal structure function [28, 50]. 
In particular, in Fig. 1 we compare the distinct predictions with the recent HI preliminary 
Fl data at Q^ = 0.75, 1.35 GeV^, as determined from the 1999 minimum bias and the 2000 
shifted vertex HI data. These points were extracted from the plots on Ref. [28]. The data 
show that Fl remains non-zero down to the lowest Q^ values measured and already distinguish 



between the different models in the low- a; region. The previous F^ data [50] at larger Q"^ are 
also presented. 

In Fig. 1, the lowest Q^ bins are shown, in particular the new Q^ = 0.75 and Q^ = 1.35 GeV^ 
measurements. At this region, the two saturation models (GBW and CGC) give very similar 
results as a consequence they have a similar behavior in the transition to the saturation regime. 
Namely, their main differences are in the large virtualities region, where CGC depends on the 
BFKL anomalous dimension at the saturation vicinity as referred before. The data description is 
very consistent, mostly at the lowest Q^ points where the usual collinear approaches are unable 
to produce reliable results. For sake of comparison, a leading twist calculation is also presented. 
We have used Eq. (2) and considered three different choices for the parton distributions. The 
predictions from the GRV98 parameterization, which is obtained using the DGLAP evolution 
equation, and the EHKQS gluon function [49] , which contains corrections from non-linear GLR 
evolution equation, are not in agreement with experimental, even at large Q^. The EHKQS 
gluon distribution slows down the dependence on x, but it is not enough to reach either to the 
upper limit of error bars. Notice that at Q^ = 0.75 GeV^ a DGLAP approach is unable to give 
reliable results, thought a backward QCD evolution is possible. For the GRV98 case, we have 
extrapolated the Q^ = 0.85 GeV^ initial condition down to Q^ = 0.75 GeV^ once they are very 
closer. Similarly, for Q^ = 1.35 GeV^, we present the EHKQS prediction for Q"^ = 1.4 GeV^ 
which is the lowest available Q^ for this parameterization. On the other hand, the predictions 
for Fl using the MRST parameterization reasonably describe the HI data, which is directly 
associated to the behavior assumed for the parton distributions in the initial evolution scale 
Q^ = 1.0 GeV^. In this case, the behavior of the sea distribution is independent of the gluon 
one, with the input gluon distribution being valence-like {xg oc x^'^^), while the input sea 
distribution has a steep growth at small x {xS oc x~*^'^^). 

In Fig. 2 one presents the high Q^ analysis. Once again, the saturation models give a 
rasonable description of data, producing a milder x growth than the DGLAP results. Their 
overall normalization also is smaller than the DGLAP analysis, becoming closer as virtuality 
increases. When comparing the DGLAP results, one verifies that the intermediate Q^ region 
at small-a; is an adequate kinematical region to study non-linear QCD corrections to the gluon 
distribution. There, very precise measuments of Fl or large statistics could constraint the size 
of those corrections. 

Finally, in Fig. 3 we present a comparison among the contributions of the different twists. 
In order to do this, we show separately the result from the summation of the different twists 
for the analitically simple expression in Eqs. (5-6) for the saturation model. For comparison 
we also present the GBW prediction, which represents the full higher-twist ressumation. In 
order to compare the present analysis with a leading order and leading twist calculation, we 
present in addition the results of Eq. (2) considering the GRV98 gluon distribution. We can 
see that the twist-2 contribution from the saturation model is in complete numerical agreement 
with the DGLAP result, showing a consistent reproduction of leading twist contribution in 
its twist resummation. We have that the twist-4 and twist-6 terms gives important contribu- 
tions, strongly modifying the magnitude the longitudinal structure function and its x-behavior. 
Moreover, the x value where the curves become distinct is Q^-dependent, as expected from 
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Figure 2: The results for Fl{x, Q"^) as a function of x at fixed large Q^ values. Same notation 
of previous figure. 

the energy dependence present in the saturation scale and consequently in the variable ^. In 
particular, for Q^ = 4.2 GeV^ we have the sum of the twist-2 and twist-4 terms reasonable 
reproduces the GBW prediction for x > 10~^, while for Q^ = 2.2 GeV^ this approximation 
is only valid for x > 10~^. For smaller values of Q^ only the full ressumation gives a good 
description of the experimental data, which demonstrate a twist summation term by term (i.e., 
summing the first contributing terms) would be either incomplete in that kinematical regime. 

As a summary, we have analyzed the longitudinal structure function at low Q^ and small-x, 
which directly depends on the gluon distribution function, within the saturation approach. In 
particular, we have show the saturation models (GBW and CGC) describe consistently the 
recent low Q^ HI data, even at very low virtualities. Moreover, we have presented the higher 
twist contributions, using the simple analytical expressions provided by the GBW model. It 
is verified that they play an important role at small-x for the virtualities considered here. 
Concerning the leading twist DGLAP analysis, which we have considered for sake of comparison, 
it is shown that non-linear GLR corrections to the gluon distribution function are not enough 
to bring the numerical analysis to the recent experimental results. However, we have shown 
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Figure 3: The results for Fl{x, Q^) as a function of x at fixed large Q^ values for the saturation 
model (GBW) and the sum of its twist contributions. The leading twist (twist-2) and the first 
2 higher twist (twist- J^ and twist-6) correction are presented. The leading twist DGLAP results 
(using GRV98 gluon pdf) is shown for comparison. 

the intermediate Q^ region should be an important kinematical region where the size of such 
corrections could be investigated. In general lines, F^ is an outstanding observable testing both 
parton saturation and twist resummation. Therefore, more precise data and/or more statistics 
is increasingly desirable at low Q^ and small-x. Significant further progress in Fl measurements 
at HERA can only be made by reducing the proton beam energy. A run with reduced proton 
beam energies is planned for the next few years. 
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